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Abstract. The next-to- leading corrections to the jet vertex which is rel- 
evant for the Mueller-Navelet jets production in hadronic collisions and 
for the forward jet cross section in lepton-hadron collisions are presented 
in the context of a kt factorizazion formula which resums the leading and 
next-to-leading logarithms of the energy. 



1. Introduction 

Recently, in the study of QCD in the Regge limit, a novel element has 
been defined and computed at the next-to-leading order (NLO). It is the 
jet vertex [1], which represents one of the building blocks in the production 
of Mueller-Navelet jets at hadron hadron colliders and of forward jets [2] 
in deep inelastic electron proton scattering. Such processes should provide 
a kinematical environment for which the BFKL Pomeron [3] QCD analysis 
could apply, provided that the transverse energy of the jet fixes a pertur- 
bative scale and the large energy yields a large rapidity interval. 

In a strong Regge regime important contributions, or even dominant, 
come, in the perturbative language, from diagrams beyond NLO and NNLO 
at fixed order in as- This is the main reason for considering a resumma- 
tion of the leading and next-to-leading logarithmic contributions as com- 
puted in the BFKL Pomeron framework. The lackness of unitarity, if the 
related corrections are not taken into account, forces one to consider an 
upper bound on the energy to suppress them. It is already known that 
the leading logarithmic (LL) analysis is not accurate enough [4], being the 
kinematics selected by experimental cuts far from any asymptotic regime. 
Moreover at this level of accuracy there is a maximal dependence in the dif- 
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ferent scales involved (renormalization, collinear factorization and energy 
scales). For the Mueller-Navelet jet production process the only element 
still not known at the NLO level was the "impact factor", which describes 
the hadron emitting one inclusive jet when interacting with the reggeized 
gluon which belongs to the BFKL ladder, accurate up to NLL [5]. The jet 
vertex, now computed, is the building block of this interaction. For the so 
called forward jet production in DIS the extra ingredient necessary is the 
photon impact factor, whose calculation is currently in progress [6, 7]. Let 
us also remind that NLL BFKL approach has recently gained more theo- 
retical solidity since the bootstrap condition in its strong form, which is the 
one necessary for the self-consistency of the assumption of Reggeized form 
of the production amplitudes, has been stated and formally proved [8]. This 
relation is a very remarkable property of QCD in the high energy limit. 




t 



Figure 1. High energy process with jet production. H is the incoming hadron providing 
a parton a (gluon g/quark q) with distribution density /a which scatters with the parton 
b with production of a jet J in the forward direction (w.r.t H) and i is the generic label 
for outgoing particles. 



A theoretical challenge, interesting by itself and appearing in the cal- 
culation, is related to the special kinematics. The processes to be analyzed 
is illustrated in Fig.l: the lower parton emitted from the hadron H scat- 
ters with the upper parton q and produces the jet J. Because of the large 
transverse momentum of the jet the parton is hard and the collinear fac- 
torization allows for a partonic scale dependence described by the DGLAP 
evolution equations [9]. Above the jet, on the other hand, the kinematics 
chosen requires a large rapidity gap between the jet and the outgoing par- 
ton q: such a situation is described by BFKL dynamics. Therefore the jet 
vertex lies at the interface between DGLAP and BFKL dynamics, a situ- 
ation which appears for the first time in a non trivial way. As an essential 
result of our analysis we find that it is possible to separate, inside the jet 
vertex, the collinear infrared divergences that go into the parton evolution 
of the incoming gluon/quark from the high energy gluon radiation inside 
the rapidity gap which belongs to the first rung of the LO BFKL ladder. 
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2. Jet Vertex: definition and outline of the derivation 

Let us consider, for the process ihustrated in Fig. 1, the kinematic variables 

PH = (7^2,0,0^ , p^ = xpH , Pb= (o,^/7/2,0^ 
Pi = Ei(ey^/V2,e-y^/V2,4>i) , s:={pH+Pbf. (1) 
We study the partonic subprocess a + b ^ X + jet in the high energy limit 
Aqcd < ~ -t (fixed) < s ^ oo. (2) 

Starting from the parton model, we assume the physical cross section to be 
given by the corresponding partonic cross section da (computable in per- 
turbation theory) convoluted with the parton distribution densities (PDF) 
/a of the partons a inside the hadron H. A jet distribution Sj, with the 
usual infrared safe behaviour, selects the final states contributing to the 
one jet inclusive cross section that we are considering. We choose, as jet 
variables, rapidity, transverse energy and azimuthal angle and the one jet 
inclusive cross section initiated by quarks and gluons in hadron H is there- 
fore written as 



dcr da^H 



J2 y'dxdaba(x)Sj(a;)/W(x) . (3) 



dJ dyjdEjd(/)j ^^^^ 

At the lowest order the jet cross section, dominated by a t-channel gluon 
exchange, can be written as [1] 



= E / / dfe 4°^(fc)FW(fc,x)/f (x) (4) 



da(0) _ 

a=q,g • 

where Va^'^\k,x) = h^^\k)sf\k , x) is the jet vertex induced by parton a, 

h^^ (fc) is the partonic impact factor (generally expressed in D = 4 + 2£: 

dimensions for regularization pourposes at the NLO) and fa^\x) is the 
parton distribution density (PDF). The jet distribution is in this case trivial 

5f (fc,x) := 5f (m,P2;Pg,Pq) = 6{1- xj/x)E]+^'S{k - kj) (5) 

with xj := Eje^-' / y/s. 

In the NLO approximation virtual and real corrections enter in the 
calculation of the partonic cross section data- The three partons, produced 
in the real contributions, in the upper, and lower rapidity region are denoted 
by 2 and 1 respectively, while the third, which can be emitted everywhere 
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in rapidity, by 3. Moreover we shall call k = — P2 and k' = pi — Pa and 
q = k — k' . A Sudakov parametrization is chosen 

k = ~wpg + wpq + k± , k± = (0, 0, k) (6) 
k' = -zpg + zpq + k'^ , k'^ = (0, 0, k') (7) 

with the typical inequalities i(;~z<Cu)~z<Cl valid in the multi 
Regge kinematics (MRK). Dimensional regularization is used, as usual, to 
trace the infrared and ultraviolet divergences. On taking into account the 

(3) 

infrared (IR) properties of the jet distribution Sj , the following structure 
is matched exactly up to NLO (i.e. a^) [1] 



da 
dJ 



= J2 [ dx f dkdk' hb{k)G{xs,k,k')Va{k',x)fa{x) , (8) 



where 



h = + a,/i(i) + • • • , 
y = y(o)+a^y(i) + ... , 

/ = /(°) + as/(^) + • • • , 

G{xs, k, k') := 5{k - k') + a^K^^^ (fe, k') log — + • • • (9) 

So 

The partonic impact factor correction in forward direction ^^^^ is well 
known [10], the PDF's /a are the standard ones satisfying DGLAP evo- 
lution equations 

as/a«(x, Ml) := g ^ (f ) ' E ^ab /r (10) 

with jip the collinear factorization scale, and the BFKL Green's function G 
is defined by the LO BFKL kernel K'^^^ . The new element is the correction 
to the jet vertex V^^^ whose expressions can be found in [1]. 

Let us note that the previous form in (8) is clearly suggested by the 
structure of the leading logarithmic part, but at this level it is not a trivial 
ansatz. Only after a proper treatment of all the IR singular terms one will 
be able to extract and define the jet verteces initiated by quarks and gluons 
in the hadron. 

The NLO virtual corrections in the high energy limit are relatively sim- 
ple and they can be summarized in the following expression: 



dc7(™t) 

dJ 



a,Y.jdxjdk h^^\k) 



2a; W (fe) log ^ + nb(fc) + Ha (fe) 
k 



xs 



xyW(/c,x)/io)(x), (11) 
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where cu^^\k) is the 1-loop reggeized gluon trajectory and the Ili{k) are 
related to the part constant in energy, after having performed the ultraviolet 
(UV) e-pole subtraction occuring in the renormalization of the coupling. 

The NLO real corrections, even in the high energy limit, are much more 
involved since the interplay of different IR soft and collinear singularities 
follows a more complicated pattern. 

One starts from the partonic differential cross sections dabs{x), appear- 
ing in (3), which have been also computed in the high energy regime [10], 
since in such a case all terms suppressed by powers of s may be neglected. 
The form of the partonic differential cross section turns out to be quite 
simple when restricted to one of the two halves of the phase space, which 
arc obtained introducing the rapidity y' = y + log ^ "measured" in the 
partonic center of mass frame, and splitting the phase space into two semi- 
spaces defined by y'l, y'^ > {lower half) and y'^ < < y[ {upper half). 

In the quark initiated case one obtains for both regions a reasonable 
simple expression, while for the gluon initiated part it is convenient to 
consider separately the qq and gg final state cases. 

The "upper half region" , with y'^ < 0, leads to the rederivation of the 
partonic b-impact factor, already well known at NLL. The "lower half re- 
gion" , wherein y'^ > 0, which corresponds to z > Zcut '■ = is precisely 
the one from which we are able to extract the jet verteces. 

A very important ingredient, as previosly mentioned, is given by the 
jet definition Sj^\ As usual it selects from a generic n-particle final state 
the configurations contributing to our one jet inclusive observable and it 
must be IR safe. The last requirement simply corresponds to the fact that 
emission of a soft particle cannot be distinguished from the analogous state 
without soft emission. Furthermore, collinear emissions of partons cannot 
be distinguished from the corresponding state where the collinear partons 
are replaced by a single parton carrying the sum of their quantum numbers. 

The extraction of all the IR singularities from the real contributions 
and the matching with the virtual ones is done on employing the standard 
subtraction method. This consists in approximating the amplitudes in any 
singular region in order to extract the exact analytic singular behavior, and 
in leaving the remaining finite part in a form suitable for numerical integra- 
tion, since normally a full analytical computation is tecnically impossible. 
The analytic singular terms of both virtual and real contributions are there- 
fore combined. Moreover, with respect to standard NLO jet analysis, one 
has also to deal with the leading logarithmic terms which are entangled to 
the other IR singular terms. For the process considered it can be shown [1] 
that finally one is left with the collinear singularities associated to the im- 
pact factor of the parton b, to the definition of the distribution function 
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for the parton a in H and to the BFKL kernel term. The finite remaning 
parts correspond to the impact factors and jet verteces. 

Wc recall that the jet distribution functions become essential in disen- 
tangling the collinear singularities, the soft singularities, and the leading 
log s pieces when the particle 3 is a gluon. The same basic mechanism may 
be used for both the quark and the gluon initiated cases. 

— When the outgoing parton 1 is in the collinear region of the incoming 
parton a, i.e., yi — > oo, it cannot enter the jet; only gluon 3 can thus be 
the jet, ys is fixed and no logarithm of the energy can arise due to the 
lack of evolution in the gluon rapidity. No other singular configuration 
is found when J = {3}. 

— In the composite jet configuration, i.e., J = {1,3}, the gluon rapidity 
is bounded within a small range of values, and also in this case no log s 
can arise. There could be a singularity for vanishing gluon 3 momen- 
tum: even if the 1 || 3 collinear singularity is absent, wc have seen that, 
at very low z, a soft singular integrand arises. However, the divergence 
is prevented by the jet cone boundary, which causes a shrinkage of the 
domain of integration ~ for z — > and thus compensates the growth 
of the integrand. 

— The jet configuration J = {1} corresponds to the situation wherein 
gluon 3 spans the whole phase space, apart, of course, from the jet re- 
gion itself. The LL term arises from gluon configurations in the central 
region. Therefore, it is crucial to understand to what extent the differ- 
ential cross section provides a leading contribution. It turns out that 
the coherence of QCD radiation suppresses the emission probability 
for gluon 3 rapidity ya being larger than the rapidity yi of the parton 
1, namely an angular ordering prescription holds. This will provide the 
final form of the leading term, i.e., the appropriate scale of the energy 
and, as a consequence, a finite and definite expression for the one-loop 
jet vertex correction. 

It is therefore clear that the energy scale sq associated to the BFKL 
rapidity evolution plays a crucial role. The calculations show a natural 
choice, due to angular ordered preferred gluon emission and the presence of 
the jet defining distribution, which is also crucial to obtain the full collinear 
singularities which factorize into the PDF's. In particular, on considering 
the term giving the real part of the LL contribution, one can see that outside 
the angular ordered region 




l-z 



El 



<^ 6'3 > 01 <^ ys < yi 



(12) 



there is a small contribution to the cross section. This ordering sets the 
natural energy scale which, for example, leads to the expression of the jet 
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vertex for the case so{k,k') := (|fc'| + + \q\). A mild modification 

of such a scale can be performed without introducing extra singularities, 
contrary to what happens for a generic choice. In any case the use of a 
different scale requires the introduction of modifing terms. For the useful 
symmetric Regge type energy scale sr = \k\\k'\, one has 



G{xs,k,k') = {l + a^HL) 



1 + asK(°) log 



xs 



k\\k'\ 



(l + asHn), (13) 



where HL{k,k') = -K(-^\k,k') log i\k\ + \q\)/\k\) =HR{k',k). 

We do not present the final vertex expressions here, they can be found 
in [1]. 

3. Remarks and final formulas 

In order to extend the results of the two loop calculations, which have al- 
lowed to obtain the full NLL partonic cross sections, to our inclusive jet pro- 
duction case, it is necessary to check the relation between the partonic im- 
pact factors and the jet verteces. More precisely one may define the one loop 
correction to the "jet impact factor" as [V * = F^^) * /(°) + * /(^). 
This latter object and the corrections to the partonic impact factor h^^^ 
are different in one important aspects, i.e. their collinear singularities. Es- 
sentially the IR singularity associated to the LL real term is absent (sub- 
tracted) in the partonic impact factors, while this does not happen in the 
"jet impact factors" due to the constraining presence of the jet. In the "jet 
impact factors" are therefore hidden the full collinear singularities which 
should be factorized to obtain the right behavior of the PDF's satisfying 
the DGLAP evolution equations. 

Apart from this important difference there is a full correspondence, as 
shown in [1], between these two objects. In particular the same angular 
ordering mechanism, which sets a reference for the energy scale, is present. 
Thus any tuning of the BFKL kernel energy scale is done precisely in the 
same way for both impact factors and jet vertices and, thanks to the Regge 
factorization property valid up to NLL in QCD, one is allowed to use the 
known form of the NLL BFKL Green's function setting the scale to sr, as 
previously discussed. 

From a two loop analysis, up to af order, one expects to find the con- 
tributions 

J_^^'^ = i log2 ^ /,(o) ^(0)^(0)^(0) ^(0) + 1 ^ [/,(i)i^(o)F(o) /(O) 
a| d J 2 Sr srI 



8 



G. P. VACCA 



where the first term on the RHS is LL while the terms collected in the 
square brackets are NLL. 

Consequently, to obtain the jet cross section with accuracy up to NLL 
terms, one has to consider the NLL BFKL kernel 

K = a^K^^^ + alK^^^ (15) 

which has been computed with the Regge type scale sr = |fci||A;2l and leads 
to a Green's function, which resums all terms up to NLL in (8), 

G{xsMM) = j ^^(^^^ {ki\{l + a,HL)[u-K]-\l + a,HR)\k2). 

(16) 

The formula for the Mueller-Navelet jets [1] can be easily derived sym- 
metrizing the formula (8) for the two jet case. More precisely one obtains 



dV 

dJidJ2 , 



^ j dxi(lx2 j dfcidfc2 /a(2;i)V;(fci,xi) 

a,b 

G{xiX2S, ki, k2)Vb{k2, X2)fb{x2) (17) 



where a, b = q, g, and the subscripts 1 and 2 refer to jet 1 and 2 in hadron 
1 and 2, resp. Once the definition of the jets is given, all the elements are 
known and a definitive computation may be carried on. 

Regarding the Green's function, two more remarks should be made. 
Since the high energy asymptotic regime is far, a few iterations of the 
kernel instead of a full resummation, as given in eq. (16), may be of some 
interest in a phenomenological application. The known large corrections of 
the BFKL kernel are related to the presence of large logs. This problem 
has been solved using optimal renormalization scale methods [11] or by 
considering an improved version of the kernel, equivalent at NLL level, but 
including correcting higher order terms, following a coUinear prescription 
[12]. 



4. Conclusions 

A brief review of the definitions and calculations at the NLO accuracy of the 
jet vertex relevant for the Mueller-Navelet jets at hadron-hadron colliders 
and for forward jets in deep inelastic electron-proton scattering has been 
given. This object, defined for high energies, due to specific kinematics, 
lies at the interface between DGLAP and BFKL dynamics. Being known 
all the ingredients, it is now possibile to perform a numerical analysis of 
the production of Mueller-Navelet jets at NLL. For the forward jet pro- 
cesses further theoretical studies on the NLL photon impact factors are 
still required. 
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